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ABSTRACT 

We present deep, sub-horizontal branch, multi-colour photometry of the Andromeda II 
dwarf spheroidal (And II dSph) taken with the Subaru Suprime-Cam wide field cam- 
era. We identify a red clump population in this galaxy, the first time this feature has 
been detected in a M31 dSph, which are normally characterized as having no signifi- 
cant intermediate age populations. We construct radial profiles for the various stellar 
populations and show that the horizontal branch has a nearly constant density spatial 
distribution out to large radius, whereas the reddest red giant branch stars are cen- 
trally concentrated in an exponential profile. We argue that these populations trace 
two distinct structural components in And II, and show that this assumption provides 
a good match to the overall radial profile of this galaxy. The extended component 
dominates the stellar populations at large radius, whereas the exponential compo- 
nent dominates the inner few arcminutes. By examining colour-magnitude diagrams 
in these regions, we show that the two components have very different stellar popula- 
tions; the exponential component has an average age of ~ 7— lOGyrs old, is relatively 
metal-rich ([Fe/H] ~ —1) but with a significant tail to low metallicities, and possesses 
a red clump. The extended component, on the other hand, is ancient (~ 13Gyrs), 
metal-poor ([Fe/H] ~ —1.5) with a narrower dispersion <r [Fe/H] — 0.28, and has a well 
developed blue horizontal branch. The extended component contains approximately 
three-quarters of the light of And II and its unusual density profile is unique in Lo- 
cal Group dwarf galaxies. This suggests that its formation and/or evolution may have 
been quite different to other dwarf galaxies. The obvious chemo-dynamical complexity 
of And II lends further support to the accumulating body of evidence which shows 
that the evolutionary histories of faint dSph galaxies can be every bit as complicated 
as their brighter and more massive counterparts. 

Key words: galaxies: dwarf — galaxies: individual (Andromeda II) — Local Group 
— galaxies: stellar content — galaxies: structure 



1 INTRODUCTION 

The classic view of dwarf spheroidal (dSph) galaxies consist- 
ing of a single, old, stellar population has changed drastically 
in recent years, as deeper and more detailed observations of 
the Milky Way (MW) satellites have been conducted. The 
shift has been so dramatic that only a minority of the MW 
dSphs are now suspected of being composed of a single stel- 
lar population, with most demonstratin g multiple epochs of 
star formation (for recent reviews, see iGrebell 119971 : iMateol 



* Based on data collected at Subaru Telescope, which is operated 
by the National Astronomical Observatory of Japan 



1 19981: Skillman] | 2005l . and references therein) . As emphasised 
m IGrebell (|1997l ). no dwarf galaxies in the Local Group ap- 



pear to share the same star formation history (SFH). 

Most recently, various groups are independently show- 
ing that the evolutionary complexities revealed in the SFHs 
of some dw arf galaxies also exten d to their global structural 
properties. lHarbeck et al.l (|200ll ) investigated the presence 
of population gradients in the stellar populations of a se- 
lection of Local Group dwarf galaxies, by comparing the 
ratio of blue-to-red horizontal branch (HB) stars, and blue- 
to-red red giant branch (RGB) stars. These colour differ- 
ences should reflect changes in the ages and/or metallicities 
of the stellar populations. They found evidence for gradi- 
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ents in six out of the nine syste ms they studied. R ecent 
spec troscopic studies of Sc ulptor l|Tolstov et al.l 2004 ), For- 
nax ()Battaglia et al.ll20olih and Canes Venatici ( Ibata et al.l 
2006) have shown that these galaxies possess spatially and 
kinematically distinct stellar populations. The radial gradi- 
ents in these systems therefore reflect the changing contri- 
butions of the d istinct stellar comp onents as a function of 
radius. However. iKoch et al.l (|200Tbl Jah do not detect any ra- 
dial gradients or kinematically distinct populat ion in their 
spectr oscopic studies of Leo I and Leo II, and IKoch et al.l 
(2006) do not spectroscopically detect a metallicity gra- 
dient in Carina desp ite a strong age gradient implied by 
lHarbeck et al.l l|200ll ). At this time, it is therefore unclear 
how common all these features are, particularly regarding 
kinematically distinct populations. 

Less is known about the global stellar populations of 
the M31 dSphs than for the M W dwarfs. Howe v er, se veral 
of the dSph galaxies studied by lHarbeck et al.l (120011 ) be- 
long to the M31 su bgroup. Hubble S pace Telescope (H ST) 
WFPC2 imaging bv lPa Costa et all (|l996l . l2000l , 12002ft and 
related programs were us ed to explore the pre sence of gra- 
dients in these systems. lHarbeck et al.l l|200ll ) found that 
only Andromeda (hereafter And) I and VI in this subgroup 
show clear evidence for gradients. This agrees with the anal- 
yses by Da Costa et al. on the same data, where no radial 
gra dient was detected over the fields sampling And II and 
III (|Da Costa et al.ll2000l,l2002l). where as a gradient was ob- 
serv ed in And I (|Pa Costa et al.lll996l ) . 

iMcConnachie fc Irwir] (1200681 ) (hereafter MI6) have 
recently presented global wide-field photometry of An- 
dromeda I, II, III, V, VI and VII taken with the Isaac New- 
ton Telescope Wide Field Camera (INT WFC), sampling 
the top few magnitudes of the RGBs of these systems. By 
analysing the global stellar structure of these galaxies, they 
find that And II shows a factor of two excess of stars in the 
central regions, above that obtained by a simple extrapo- 
lation of the outer surface brightness profile to small radii. 
MI6 suggest that this is evidence showing And II consists 
of two structural components with distinct spatial distri- 
butions, similar to Sculptor and Fornax. If the stellar pop- 
ulations of these components are different, then we might 
expect radial gradients to be detected in this galaxy as the 
relative contribution of stars from each component changes. 
This gradient would pr obably act on a s cale la rger than the 
HST field observed by |Pa Costa et al.1 l|2000l ). otherwise it 
would likely have been detected by this study. The existence 
of multiple p opulations in And II w ould be consistent with 
the fact that |Pa Costa et alj |2000l ) were unable to match 
the HB morphology of this galaxy using a stellar population 
with a single age, and implied that populations of multiple 
ages must be present. 

Multi-colour photometry which reaches the HB level 
can be a ke y discriminant of str uctural variations in a dwarf 
galaxy (eg. lHarbeck et aT1l200ll ). and is a particularly pow- 
erful tool when the data covers the entire pro jected area of 
the dwarf so that a global vi ew is obtained (eg. Tolstoy et al.l 
|2004 iBattaglia et al.ll2006l ). In addition, data of this type 
provides much stronger constraints on the ages, metallici- 
ties and SFHs of the stellar populations in the galaxy than 
can be achieved by analysis of the RGB alone. In this pa- 
per, we present deep, global VI photometry of And II which 
reaches to below the HB level, obtained as part of a survey of 



Local Group dwarf galaxies with the Subaru Suprime-Cam 
wide field camera. We use these data to probe the structures, 
ages and metallicities of the stellar populations of this dwarf 
galaxy over its entire spatial exten t, to a depth equ i valent to 
the earlier HST-WFPC2 study bv lPa Costa et al.l l|2000l ). 

This paper is organised as follows: in Section 2, we in- 
troduce the survey, and discuss the targets, observations and 
data analysis procedure. In Section 3, we present colour 
- magnitude diagrams (CMPs) for And II, and examine 
the spatial properties, ages and metallicities of the domi- 
nant stellar populations. In Section 4, we discuss the struc- 
ture and evolution of And II in light of these results. 
Section 5 summarises. We assume a distance modulus to 
And II of (m-M) = 24.07 ± 0.06 (d = 652 ± 18kpc; 
IMcConnachie et al.1 12004 120051 ). The averag e extinction in 
the d irection of And II is E (B - V) = 0.063 (|Schlegel et all 
GUI). 



2 THE SUBARU SUPRIME-CAM SURVEY 

Inevitably, the majority of our detailed information on the 
stellar content and evolution of dSphs comes from ob- 
servations of the MW subgroup, although there are sev- 
eral notable exceptions. Environmental effects, such as ram 
pressure stripping and tidal effects, are thought to play 
an important role i n the evolution of d warf galaxies (eg. 
lEinasto et all 1 19741 ; iMaver et all l2001al lbl). and so dwarfs 
which have evolved in different environments could possess 
very different properties. For example, MI6 have shown that 
the dSph satellites of M31 are a factor of 2 — 3 times more 
extended than for the MW population, showing that, in at 
least one respect, results derived from the MW population 
do no t necessarily hold for o t her Loc al Group dSphs. Sim- 
ilarly, IMcConnachie fc Irwinl (2006b) show that the radial 
distribution of all of M31's dwarf satellites is nearly twice as 
extended as that of the MW population. It is clearly of con- 
siderable importance to obtain a complete census of nearby 
dwarf galaxy properties, but it is only relatively recently 
that the instrumentation has existed to probe the global 
stellar content of the more distant dwarf galaxies of the Lo- 
cal Group. 

As a follow- up to our INT WFC survey of L ocal Group 
dwarf galaxies (|McConnachie et all 120041 . 120051 . MI6), we 
have undertaken a deeper, wide-field, multi-colour photo- 
metric survey of Local Group dwarf galaxies using the Sub- 
aru Suprime-Cam wide-field camera, with the aim of obtain- 
ing photometry which reaches below the HB components in 
each of our targets over their entire spatial extent. Pata of 
this type are generally lacking for the galaxies of the Local 
Group which are not satellites of the MW, but hold a large 
amount of information relating to their global structural 
properties and SFHs. Subaru Suprime-Cam is ideally suited 
to this task; it is able to reach the HB of Local Group dwarf 
galaxies located beyond 500 kpc in a reasonable amount of 
time and, importantly, the large (34arcmins x27arcmins) 
field of view is sufficient to cover the entire dwarf galaxy at 
this distance. 
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Figure 1. Top panels: the V— (left) and I c — band (right) Suprime-Cam fields centered on Andromeda II (34arcmins x27 arcmins). 
North is to the top, and east is to the left. The fields are affected by the presence of scattered light from bright stars either in, or just 
off, the field of view. The scattered light patterns are complex, but are generally smoothly varying. Bottom panels: same as top, where 
the scattered light has been removed using a non-linear equivalent of unsharp masking prior to stacking. Low level artifacts and negative 
halos around bright stars still remain, but this has a negligible effect of the photometry of the vast majority of objects. 
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Figure 2. Photometric errors as a function of magnitude for 
the V (top panel) and / (bottom panel) data. The data start to 
become incomplete at a signal-to noise ratio of 10, corresponding 
to V = I ~ 25.5. 



2.1 The survey 

During the nights of 3 rd - 5 th August 2005, we obtained 
Johnson - Cousins V and J c — band imaging of the M31 
dSph galaxies Andromeda I, II, III, V, VI and VII using 
Subaru Suprime-Cam (P.I. N. Arimoto). Conditions were 
uniformly excellent, being photometric throughout and with 
typical seeing of 0.5 arcsecs. At the start of the night, we also 
obtained some multi-colour imaging of the distant transi- 
tion dwarf galaxy DDO210 and the isolated dSph in Cetus. 
The observations and resu lts for DDO210 are presented in 
iMcConnachie et all l|200fj| ). 

For And II, we exposed for a total of 2200 seconds 
in V and 4800 seconds in I c , split as 5 x 440 seconds and 
20 x 240 seconds dithered sub-exposures respectively. The 
telescope was typically offset ~ 20 arcsecs between sub- 
exposures. The exposure times were designed to reach below 
the horizontal branch to an equivalent depth in the V and 
I bands. The final stacked images have sub-arcsec seeing 
over the whole array, averaging 0.51 arcsec for the I c -band 
and 0.62 arcsec for the V-band. During stacking, the single 
frame object catalogues were used to improve the positional 
match with respect to the chosen reference image catalogue 
by applying an additional 6 constant linear solution after the 
application of the differential WCS. This improved the WCS 
solution to approximately one tenth of a pixel precision. 

Data were processed using a general p urpose pipeline 
for p rocessing wide-field optical CCD data (jlrwin fc Lewis! 
2001). Images were debiased and trimmed, and then flat- 
fielded and gain-corrected to a common internal system us- 
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ing clipped median stacks of nightly twilight flats. In ad- 
dition, the J c — band images, which suffer from an additive 
fringing component, were also corrected using a fringe frame 
computed from the entire series of I c -band exposures taken 
during the 3 nights. The top panels of Figure [TJshows the full 
34 arcmins x 27 arcmins field for our reduced V and I c — band 
Suprime-Cam images of And II. North is at the top, and east 
is to the left. 

The top panels of Figure Q] reveal the presence of scat- 
tered light which affects the final stacked images in an un- 
usual way. The light appears to be scattered from bright 
stars either in, or just off, the field-of-view. This leads to 
complex, but generally smoothly varying, background light 
patterns. However, when the component images are stacked 
in a conventional way the induced scattered light patterns 
can become disjoint, even after allowing for overall changes 
in the background during the stacking process. To compen- 
sate for this we were forced to remove the smoothly varying 
background components prior to stacking using a non-linear 
equivalent of unsharp masking. Although this worked well 
at removing the disjoint patterns, some of the more rapidly 
varying spatial components are still visible as low level ar- 
tifacts (eg. the horizontal striations below the center). The 
other drawback of unsharp masking is the low level negative 
halo induced around bright stars; however, this has a negligi- 
ble effect on the photometry of the overwhelming majority 
of individual objects. The bottom two panels of Figure [T] 
show the V and I c — band Suprime-Cam images after this 
correction has been applied. 

For each image frame an object catalogue was generated 
using the ob ject detection and parameterisation procedure 
discussed in Irwin et al.l l|2004l ). Astrometric calibration of 
the individual frames was based on a simple Zenithal polyno- 
mial model derived from linear fits between catalogue pixel- 
based coordinates and standard astrometric stars derived 
from on-line APM plate catalogues. The astrometric solu- 
tion was then used to register the frames prior to creating 
a deep stacked image in each passband. Object catalogues 
were then created from these stacked images and objects 
were morphologically classified as stellar or non-stellar (or 
noise-like). The detected objects in each passband were then 
merged by positional coincidence (within 1 arcsec) to form 
a V, I c combined catalogue. 

We cross-correlated the Suprime-Cam photometry with 
our earlier multi-colour INT WFC photometry of And II, for 
which we know the colour transformations into the Landolt 
systerrQ. This ensures our new photometry is on the same 
systems as our previous photometry. By only considering 
those objects reliably identified as stellar in all four sets of 
observations, we find 



V = V' + 0.030 (V - /) 
I = Ic - 0.088 (V - I) 



(1) 



where we now use V' to denote the original Subaru V filter. 
These transformations are identical to those der ived using 
our data for DDO210 jMcConnachie et al.l 120061 ) . Figure [2] 
shows our photometric errors as a function of magnitude for 
our l^-band (top panel) and /-band (bottom panel) data. 



http:/ /www. ast.cam.ac.uk/ wfcsur/technical/photom/colours 
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Figure 4. The mctallicity distribution function (MDF) for stars 
in the top two magnitudes of the RGB with r < 10 arcmins in An- 
dromeda II (solid histogram). The dashed histogram shows the 
MDF derived for the foreground population (r > 14 arcmins), 
scaled by area. These were created by interpolating between 
13 Gyr isochrones from VandenBerg et al. (2006) with BVRI 
colour-T e j f relations as described by VandenBerg & Clem (2003). 
The (foreground corrected) mean metallicity is [Fe/H] = —1.5 ± 
0.1 with fT[p e / H ] = 0.34 dex. 



Our errors are < 0.02 mags for V < 23.6 and / < 23.4. The 
data start to become incomplete at a signal-to-noise of 10, 
corresponding to V = I ~ 25.5. 



3 STELLAR POPULATIONS AND 
STRUCTURE 

3.1 Colour-magnitude diagrams 

The left panels of Figure [3] show the extinction-corrected 
CMDs for And II, where each star has had its V and I 
magnitude corrected for extinction b y cross-correlat i ng th e 
position of the star with the maps of ISchlegel et alj l|l998h . 
The right panels show the corresponding Hess diagrams 
with square-root scaling. Also marked on the CMDs are 
error bars which represent the 1 a photometric uncertain- 
ties in our data. These data go several magnitudes deeper 
than our earlier INT WFC imaging in both filters, to a 
depth compara b le to the earlier HST-WFPC2 imaging of 
|Pa Costa et al.l (|2000l ) . Suprime-Cam has a ~ 100 larger 
field of view than HST-WFPC2; > 35000 stars are shown 
in the CMDs in Figure [3] compare d to the ~ 2000 shown 
in Figure 5 of |Pa Costa et al.l (^POOp . A MW foreground se- 
quence is visible in Figure [3] at bright magnitudes, starting 
at (V -I) ~ 0.6 mags. All the other features in the CMD 
are produced by stellar populations intrinsic to And II. The 
lack of any obvious bright, main-sequence stars or blue-loop 
stars suggests that there has not been any recent (last few 
Gyrs) star formation. Instead, all of the prominent popula- 
tions (the RGB and core helium burning stars) are signs of 
intermediate (2 — 10 Gyrs) and/or old (> 10 Gyrs) ages. 
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Figure 3. Top panels: Extinction-corrected colour magnitude diagrams (left panels) and Hess diagrams (right panels) for Andromeda II. 
Error bars show the average 1 a uncertainties in the photometry at each magnitude level. The dashed boxes in the V-band Hess diagram 
show the colour — magnitude cuts used to define the loci of various stellar populations. Star symbols show the (< V >, < (V — I) >) 
locus of RR Lyrae stars in Andromeda II identified in Da Costa et al. (2000) and Pritzl et al. (2004). 



3.1.1 The red giant branch 

In the top left panel of Figure 3 we have divided the broad 
RGB of And II into a 'blue' and 'red' component. In general, 
blue RGB stars are likely younger and/or more metal-poor 
than redder RGB stars, although the unknown effects of the 
well known age-metallicity degeneracy leads to difficulty in 
the interpretation of RGB colour. We also note the presence 
of an overdensity of stars on the RGB at V ~ 23.9, which is 
most likely due to a RGB bump. 

iDa Costa et alJ l|2000t ) showed that And II has the 
broadest colour dispersion in its RGB of any of the Lo- 



cal Group dSph galaxies, implying trp wH] ~ 0-3 6 dex for 
a single, co-eval, ancient population. ICote et al.l l|l999bf ) 
spectroscopically measured a scatter in the metallicity of 
50 RGB stars in And II of <7 [Fc/H ] = 0.34 ± 0.1 dex, vir- 
tually identical to the photometric measurement. Figure [4] 
shows a metallicity distribution function (MDF) for r < 
10 arcmins derived from our data, where we have calculated 
the metallicity of each star in the top two magnitudes of 
the RGB by interpolating its position in colour-magnitude 
space between 13 Gyr isochrones with a range of metallici- 
ties (assuming zero a— enhancement). Our technique is stan- 
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Figure 5. Top panel: an enlargement of the V-band Hess diagram 
of Andromeda II in the region of the horizontal branch (HB). 
Star symbols show the (< V >, < (V — I) >) locus of RR Lyrae 
stars in Andromeda II identified in Da Costa et al. (2000) and 
Pritzl et al. (2004). The red HB begins immediately to the right 
of the RR Lyrae locus. However, the morphology of this feature 
appears to change significantly as a function of colour. The middle 
and bottom panels show luminosity functions in the two strips 
indicated: strip A will only sample red HB stars, whereas strip 
B will additionally sample some of the RGB and a red clump, 
if present. The morphology of the peaks in the two luminosity 
functions are very different, implying cither that the spread in 
red HB luminosity becomes much larger to the red, or that an 
additional stellar population (most likely a red clump) is present 
in B but not in A. Note that a RGB bump is also visible in these 
panels at V ~ 23.9. 



dard, and uses th e Vict oria-Retina set of isochrones from 
IVandenBerg et aljj ;2006) with BVRI c olour - T e f j relations 
as described by IVandenBerg fc Cleml (|2003h . The dashed 
histogram shows the MDF for RGB stars at r > 14arcmins, 
scaled by area, which will be dominated by foreground stars 
and which we use as a reference field. Since these stars are 
predominantly foreground dwarfs in the Milky Way halo, the 
metallicities which we calculate are physically meaningless, 
and useful only as a comparison to the main MDF. 

The mean metallicity of RGB stars under our age as- 
sumption is [Fe/H] = —1.5 ± 0.1. The dispersion in metal- 
licity is gpwH] = 0.35 ± 0.1 dex, in exc e llent a greement with 
|Pa Costa et al] (120001 ) and lCote et~aH 11199901 ). demonstrat- 
ing good consistency with these earlier studies. We note, 
however, that the mean metallicity and dispersion are sen- 
sitive to our age assumption, reflecting the age-metallicity 
degeneracy in the colour of the RGB discussed above. If the 
stellar populations of And II are younger than 13 Gyrs, or if 
a range in age is present, then the metallicity estimates we 



derive will be misleading. Later, we show that this is nearly 
certainly the case for And II. 

3.1.2 The core helium burning stars 

Horizontal branch stars (low mass, core helium burning) are 
only present in old stellar populations. The presence of blue 
HB stars and RR Lyrae variables in particular are an unam- 
biguous sign of the presence of a population which is at least 
as old as the MW globular clusters (> 10 Gyrs). The CMDs 
in Figure 3 clearly show the presence of a HB at V ~ 24.7 
which extends far to the blue. The boundary between blue 
and red HB stars is marked by the instability strip: at the 
HB level, this is marked by the mean colours of RR Lyrae 
variable stars. 

The Subaru Suprime-Cam dat a is unsuitable for the 
identification of RR Lyraes. However. |Pa Costa et all (|2000h 
identify over 70 RR Lyrae variables in their HST-WFPC 2 
dataset, which are analysed in detail in IPritzl et al] i|2004f ). 
and who give the < V > and < (B — V) > magnitudes for 
these stars. To compare the position of this RR Lyrae lo- 
cus with our V, I data requires the calculation of the colour 
transformation between (B - V) for the HST-WFPC2 data 
and (V — I) for the Suprime-Cam data. To this end, we have 
cross-correlated the HST-WFPC2 dataset with the Suprime- 
Cam data and identified common objects by positional coin- 
cidence within 1 arcsec. We apply the additional constraint 
that the photometric uncertainties in the ground-based data 
should be less than 0.05 magnitudes in each filter. We then 
determine the relation between (B — V) and (V — 7), and 
find it is well fit by the linear relation 

(V -I) = 0.837(75-1/)+ 0.175 . (2) 

The star symbols in the top right panel of Figure 3 show the 
(< V >, < (V — 7) >) locus of the RR Lyrae stars identified 
in the earlier HST-WFPC2 study. The rectangle surround- 
ing these points shows the colour-magnitude cuts we use to 
approximate the position of this locus. The large rectangle 
immediately to the blue of this defines the colour-magnitude 
cuts used to isolate the blue HB. Likewise, the red HB will 
lie immediately to the red of the RR Lyrae locus. However, 
the structure of the CMD in this region is quite complex. 

The top panel of Figure [5] shows an enlarged version 
of the V— band Hess diagram at the level of the HB. The 
RR Lyrae and blue HB are marked as before. To the red 
of the RR Lyrae locus, a significant overdensity of stars is 
present. Some of these will clearly be red HB stars; how- 
ever, the width of this overdensity increases significantly to- 
wards the red. To illustrate this, in the middle panel we 
have plotted the luminosity function of stars in a strip in 
colour space immediately redder than the RR Lyrae locus 
(A) and compared this to the luminosity function of stars in 
another strip which is even redder (B) in the bottom panel. 
Strip A is likely to only sample red HB stars; however, strip 
B additionally sample a significant number of RGB stars, 
and will also sample red clump stars, if present. This lat- 
ter stellar population are also core helium burning stars, 
but they have a significantly higher mass than their low 
mass, HB counterparts. Their presence would therefore im- 
ply the existence of an inte rmediate-age stellar population 
(see lGirardi fc SalarisllioOll for a comprehensive theoretical 
study of the red clump). 
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The peak of the luminosity functions in Figure [5] corre- 
spond to core helium burning stars in And II; however, the 
peaks clearly have vastly different dispersions and shapes. 
Therefore, either the spread in the luminosity of red HB stars 
varies dramatically as a function of colour, or the broader 
peak in B is due to the presence of an additional stellar 
population with similar luminosities and colours to the HB. 
We are unaware of a physical mechanism w hich could cause 
the fo rmer to occur. On the other hand, iDa Costa et al.l 
(2000) independently implied that And II had to possess 
an intermediate-age population (of age 6 — 9 Gyrs) in order 
to explain the overall HB morphology of And II as revealed 
by their HST-WFPC2 data. Such a population would be ex- 
pected to possess a red clump. Therefore, we conclude that 
the morphology of the region of the CMD shown in Fig- 
ure [5] is best and most naturally explained by the presence 
of a red clump population in this galaxy. This is the first 
detection of a red clump in any of the M31 dSphs, which 
are often characterised as having no significant intermediate 
age populations. Colour-magnitude cuts to define the red 
HB locus and red clump locus are shown in the top-right 
panel of Figure 3 (the red clump locus will also contain sig- 
nificant contributions from red HB and RGB stars). We will 
examine the luminosity and colours of all the core helium 
burning stars in more detail in Section 3.5. 



3.2 Radial gradients 

3.2.1 The red-to-blue horizontal branch profile 

IDa Costa et~aH l|2000r i and lHarbeck et al l |200lT l looked for 
population gradients in And II by examining the change in 
the number of blue HB stars to the total number of (blue 
+ red) HB stars over the radial extent of the HST-WFPC2 
field. Both papers concluded that no significant radial vari- 
ation in this quantity was present in their data. 

Using our colour cuts to define blue an d red HB stars, 
we ha ve calculated an analogous quantity to lDa Costa et aT] 
(2000), Ni — urhb I tibhb (the ratio of the number of red 
HB stars to the number of blue HB stars) as a function of 
radius. The radial extent of And II probed by our data is 
significantly larger than that probed by the HST-WFPC2 
data, and so should display any large-scale variations in 
this ratio if they are present. We count stars in annuli of 
width 1 arcmin centered on And II, and display our results 
in the top panel of Figure [6] JVi remains approximately flat, 
with no significant deviation from the mean value, out to 
13 arcmins from the center of And II. This implies that all 
of the obviously old stellar populations in And II have very 
similar radial distributions, and is in agreement with the 
earlier studies. 



3.2.2 The red clump to horizontal branch profile 

We now define an analogous quantity to iVi , Na = urc /ns b 
(the ratio of the number of red clump stars to the number 
of all HB stars). This will probe the relative contribution 
of low-to-high mass helium burning stars in And II as a 
function of radius. 

The profile of N2 is shown in the bottom panel of Fig- 
ure [6] In contrast to the ratio of red to blue HB stars, the 
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Figure 6. Top panel: the ratio of the number of red HB stars 
to the number of blue HB stars, as a function of radius in An- 
dromeda II. The dotted line is the mean value of this ratio aver- 
aged over all radii. Clearly, this ratio does not change significantly 
over a large radial range, implying that the blue HB and red HB 
stars posses similar radial distributions. Bottom panel: the ratio 
of the number of red clump stars to the number of (blue + red) 
HB stars. The ratio decreases significantly as a function of radius, 
demonstrating that these populations possess different spatial dis- 
tributions, such that the red clump stars are significantly more 
concentrated than the HB stars. 



ratio of red clump to HB stars declines significantly and con- 
tinually as a function of radius by nearly a factor of 6 over 
a 13 arcmin radial range. Clearly, the low and high mass 
helium burning stars in And II have very different radial 
distributions, such that the red clump stars are more cen- 
trally concentrated than the HB stars. This then implies 
that the intermediate age stellar populations in And II are 
significantly more centrally concentrated than the old stel- 
lar populations. And II clearly possesses a very strong radial 
gradient in its stellar populations. 



3.3 Spatial distributions 

Figure [7] shows maps of the spatial distribution of stellar 
sources satisfying various colour-magnitude cuts in And II, 
projected into the tangent plane of this galaxy. The top 
panels of Figure [7] show the spatial distribution of all stars 
in And II brighter than our ~ 90 % incompleteness limits. 
And II is a very large, extended spheroid of stars with a 
slight ellipticity of e ~ 0.2 (MI6). Foreground stars will dom- 
inant in the outer parts of the Suprime-Cam field, but the 
tidal extent of And II is sufficiently large (r< = 19.8 arcmins; 
MI6) that some of its stars will be present at all radii probed. 

The remaining five panels in Figure [7] correspond to 
stars in the stellar loci indicated in Figure 3. Given the tight 
colour-magnitude cuts used to isolate each population, the 
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Figure 7. The spatial distribution of stellar populations in Andromeda II. The top right panel consists of all stars brighter than 
our ~ 90 % completeness limits. The five remaining panels show the distribution of various stellar populations, defined by the colour- 
magnitude loci shown in Figure 3. The blue and red HB populations both posses diffuse and extended distributions. In contrast, the red 
RGB stars have a very different spatial distribution which is significantly more concentrated and less extended than the HB distribution. 
The blue RGB and red clump distributions arc relatively similar to each other, but are different to the other sub-populations. 



contamination from foreground stars in these panels is gen- 
erally small. Some of the irregularity of the stellar distribu- 
tions can be traced back to the presence of saturated stars 
in the field, which can cause holes to appear in the spatial 



distributions; for example, near the center of the red and 
blue HB distributions. 

The red and blue HB both have extended, diffuse stellar 
distributions which do not obviously increase in density at 
small radius, as is common for spheroidal profiles. In com- 
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parison, the red clump and blue RGB populations are ap- 
proximately as extended as the HB distributions, but are 
much more centrally concentrated. This naturally explains 
the gradient observed in Section 3.2.2. The red RGB, on the 
other hand, is far less extended than all the other stellar 
distributions, and is concentrated in the central regions of 
And II. The contrast between the HB distribution and the 
red RGB distribution is striking. 

3.4 Radial profiles 

We robustly quantify the spatial distribution of the stellar 
populations shown in Figure [7] by constructing radial pro- 
files for each population. We count stars in elliptical annuli 
with the average ellipticity and position angle of And II 
(e = 0.2,6* = 34°; MI6). While the maps in Figure [7] show 
that these quantities change for And II depending upon the 
stellar population examined, adopting annuli of fixed shape 
and orientation for all the populations allows for a more ro- 
bust and meaningful comparison between the resulting pro- 
files. 

The radial profiles for the individual stellar populations 
are shown in the left and center panels of Figure [8] in log-log 
space. The HB profile consists of the blue HB, RR Lyrae and 
red HB stellar loci, since these have been shown to possess 
similar distributions. The radial profile for all stars above 
the ~ 90 % completeness limit of our data is shown in the 
right panel of Figure [S] 

A foreground correction has been applied to the overall 
profile by comparison with the (foreground-corrected) radial 
profile in MI6. We have assumed that the core and tidal 
radii of the best-fitting King model derived by MI6 should 
also fit the outer regions of our Suprime-Cam data. Our 
free parameters in this fit are a vertical scale factor and a 
constant representing the foreground contamination in the 
Suprime-Cam data. The best-fit value of the latter can then 
be subtracted from the Suprime-Cam profile. Although no 
foreground correction has been applied to the profiles for 
the individual populations, the tight colour-magnitude cuts 
which define these populations ensure that the foreground 
contamination is minimal. 

Since the radial profiles in the left and center panels 
of Figure [5] show significant differences, we choose the flex- 
ible Sersic profile to parameterise them all. The projected 
density of this distribution function is given by 

/(r)=/ exp[-(r/r o ) 1/n ] . (3) 

n is a shape parameter (n = 1 corresponds to an exponential 
profile and n — 4 corresponds to a de Vaucouleurs profile) 
and r is a scale-radius. To is a normalisation constant and 
is equivalent to the central surface brightness. The best-fit 
Sersic profiles to the radial profiles of the individual stellar 
populations are shown as dashed lines in the left and center 
panels of Figure[8] For each population, we find the following 
best-fit parameters: 

• Red clump: r a = 5.76 ± 0.84, n = 0.60 ± 0.02 

• Blue RGB: r a = 4.34 ± 0.47, n = 0.79 ± 0.02 

• Red RGB: r a = 2.18 ± 0.78, n = 1.03 ± 0.29 

• HB: r a = 10.00 ± 8.48, n = 0.30 ± 0.11 

The red clump and blue RGB are the most similar out of 



these 4 profiles, although even these are statistically different 
from each other. The red clump profile has a slightly more 
steeply declining profile than the blue RGB, and a slightly 
larger scale length. Both these profiles are very distinct from 
the red RGB profile, which has an essentially exponential 
(n — 1) profile with a scale length of only 2.18 ± 0.78 ar- 
cmins (417 ± 152 pc). This is much more similar to the typi- 
cal overall radial profile of dSp h galaxies in the Local Group 
(|lrwin fc Hatzidimitrioul 11995. MI6). We note with interest 
that the exponential scale length of this population is a fac- 
tor of ~ 2 larger than the exponential scale length of the 
overall radial profiles of most of the MW dSph population 
jlrwin fc Hatzidimitrioull 19951 ). 

The radial profile of the HB component is very unusual. 
The slight decline in the stellar density of this component 
at small radii is an artifact of the few saturated stars in this 
region which produce holes in the spatial maps in Figure [7] 
The scale radius of the HB profile is not well constrained, due 
to the unusual, nearly constant density shape of the profile. 
The value of the power-law index, n = 0.30±0.11 produces a 
very steeply declining profile at large radius. Thus, the HB 
population has an approximately constant density profile 
out to large radius before it essentially truncates. 

Although the radial profiles of the red RGB and HB can 
both be described by Sersic profiles, the two bear no similar- 
ity to each other at all, as even a cursory glance of Figure 8 
makes clear. A K-S test on the red RGB and HB cumulative 
radial profiles shows that the two are absolutely inconsistent 
with being drawn from the same underlying distribution at 
extremely high significance (>> 99.99%). 

Given that the red RGB has an exponential profile, 
and the HB has a constant density profile, it is very likely 
that the two must have very different dynamical proper- 
ties. For example, if they are both in equilibrium with the 
overall potential of the dSph, then they cannot both have 
the same kinematic properties (velocity dispersion profile, 
orbital anisotropy etc) and have such vastly different den- 
sity profiles. This is because they will not both satisfy Jeans 
equation, which requires a balance between the overall po- 
tential, the density distribution of the tracer population and 
the kinematic properties of the tracer population. On the 
other hand, if one or both of the populations is not in equi- 
librium, then there is no reason to expect that they will have 
similar kinematics in the first place. It is very likely there- 
fore that the red RGB population and the HB population 
belong to two kinematically distinct components in And II. 

Compelling evidence for this interpretation comes from 
the overall radial profile displayed in the right panel of Fig- 
ure[S] The dot-dashed line is the best-fit to the data obtained 
on fitting a weighted sum of the HB and red RGB profiles, 
that is the function 

F(r) = aiI H B{r) + a 2 InRGB{r) (4) 

where Ihb(t) and Irrgb{t) are the best fitting Sersic pro- 
files to the HB and red RGB populations found previously, 
and ai and a 2 are scaling constants found through a least- 
squares fit to the overall profile. If, as we postulate, the HB 
and red RGB populations trace two distinct components in 
And II, then the sum of the radial profiles of these compo- 
nents must equal the overall radial profile of And II. 

The dot-dashed line in the right panel of Figure [8] is 
a remarkably good fit to the overall profile of And II. The 
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Figure 8. Left and center panels: radial distributions in elliptical annuli of the stellar populations in Andromeda II denned by the 
colour- magnitude loci in Figure 3 (the blue HB, RR Lyrae and red HB have been grouped together). Best-fitting Sersic profiles have 
been overlaid as dashed lines. The red RGB and HB profiles are most distinct; the former is nearly exponential whereas the latter 
is approximately constant density out to ~ 10 arcmins. Right panel: the overall radial profile of Andromeda II derived using all stars 
brighter than our ~ 90% completeness limit. Remarkably, the overall profile of this galaxy is very well described by an appropriately 
weighted addition of the HB and red RGB profiles. The fit obtained on doing this is shown as a dot-dashed line, whereas the scaled HB 
and red RGB profiles are shown as dotted and dashed lines, respectively. This suggests that Andromeda II is composed of two dominant 
components; a centrally concentrated, nearly exponential component best traced by red RGB stars, and an extended, nearly constant 
density component best traced by the HB stars. 



only slight deviation of this line from the data is at large 
radius, where we are particularly sensitive to the foreground 
correction we applied (which was not measured directly from 
our data). Therefore, this deviation seems small as it is only 
of order 1 star arcmin -2 , and the fit is good over 2.5 dex in 
stellar density. 

The individual, weighted Sersic profiles of the HB and 
red RGB populations, a\InB{r) and cl2Irrgb (r) respec- 
tively, are shown in the right panel of Figure [5] as dotted 
and dashed lines respectively. Inside of 2 - 3 arcmins, the 
exponential component traced by the red RGB stars dom- 
inates, whereas outside of this radius the extended compo- 
nent traced by HB stars begins to dominate. This natu- 
rally explains the peculiar profile for And II measured by 
MI6, where a factor of two excess of stars was seen at 
r < 2 arcmins upon the extrapolation of the best-fitting den- 
sity profile at large radius to small radius: at large radius, 
the extended component dominates, whereas at small radius 
an additional component contributes to the density profile, 
and dominates over the extended component. In Section 3.5, 
we present additional evidence in the ages and metallicities 
of the stellar populations of And II which supports the pres- 
ence of two distinct components. We note that the blue RGB 
and red clump profiles are also able to be approximated by a 
weighted sum of the red RGB and HB profiles, as is required 
in this model. 

The absolute magnitude of And II obtained by MI6 



upon integrating the surface brightness profile is Mv = 

— 12.6 ± 0.2. By using the best-fit values of a\ and a,2, it 
is straightforward to calculate the expected luminosities of 
the two postulated components. We find that the extended 
component contributes ~ 74% of the light in And II, while 
the exponential component will contribute ~ 26 % of the 
light. This translates into magnitudes of My ~ —12.3 and 

— 11.2 for the extended and concentrated components, re- 
spectively. The corresponding central surface brightnesses 
are /xo,v ~ 26.1 and 24.8 mags arcsec -2 , respectively, where 
the total central surface brightness of And II has been taken 
to be 24.5 mags arcsec -2 (MI6). 

Table 1 lists the above parameters for And II, along 
with error estimates. The uncertainties in My and Lv are 
obtained by propagating the uncertainties in the absolute 
magnitude of And II, a\ and a,2. They are likely to be under- 
estimates of the true uncertainties, since the adopted profiles 
of each component also have considerable uncertainties as- 
sociated with them. Nevertheless, these values illustrate the 
approximate contribution of each component to the overall 
composition of And II. 

3.5 The ages and metallicities of Andromeda II 

The analysis in Section 3.4 suggests that different stellar 
components with different characteristic stellar populations 
dominate And II at small and large radius. To investigate 
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Figure 9. Extinction-corrected V vs (V — I) and I vs (V — I) CMDs for Andromeda II for r > 7.5 arcmins (top panels) and r < 2 arcmins 
(bottom panels). The former will predominantly sample the stellar populations of the extended component of Andromeda II, whereas 
the latter will predominantly sample the stellar populations of the centrally concentrated, exponential component. A 13 Gyr isochrone 
with [Fc/H] = —1.5, from the Victoria - Regina set of isochrones, is overlaid to aid visual comparison. The differences between the CMDs 
split in this way are striking, particularly regarding the width of the RGB and the morphology of the core helium burning stars. 



this possibility further, we show in Figure [9] the extinction 
corrected CMDs of And II split by radius. The top panels 
show the And II CMDs for stars with r < 2 arcmins. Ac- 
cording to the fits in Section 3.4, the exponential component 
dominates the extended component by a factor of approxi- 
mately 2:1 inside this radius. The bottom panels of Figure[9] 
show the And II CMDs for stars with r > 7.5 arcmins. The 
extended component is expected to dominate the concen- 
trated component by a factor of approximately 8:1 over this 
radial range. To aid visual comparison, a 13 Gyr isochrone 



with [Fe/H] = —1.5 from the Victoria - Regina set of 
isochrones, shifted to the distance modulus of And II, has 
been overlaid. The difference in the stellar populations di- 
vided by radius is striking, particularly regarding the width 
of the RGB and the morphology of the HB/red clump level. 
Each of these features contains information regarding the 
average ages and metallicities of the dominant stellar popu- 
lations at these radii. 
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Figure 10. Top panels: Extinction-corrected V and /-band luminosity functions for stars in Andromeda II with r > 7.5arcmins (which 
predominantly sample the extended component). The dashed lines represent the best-fits of a power law and Gaussian to the RGB in 
the region of the red clump/HB level. The Gaussians are centered at Vjjb = 24.76 ± 0.01 and Ijjb = 24.04 ± 0.01, with dispersions of 
ay = 0.17±0.01 and aj = 0.14±0.01. Bottom panels: Extinction-corrected V and /-band luminosity functions for stars in Andromeda II 
with r < 2arcmins (which predominantly sample the exponential component). The dashed lines represent the best-fits of a Gaussian to 
the red clump/HB level. The Gaussians are centered at VhcB = 24.73±0.02 and IhzB = 23.81±0.02, with dispersions of ay = 0.37±0.05 
and 07 = 0.34 ± 0.02. 



3.5.1 Populations at large radii 

The large majority of stars in the top panels of Figure [9] are 
expected to belong to the spatially extended stellar com- 
ponent. A blue HB population occupying the same colour- 
magnitude locus as denned earlier is clearly visible in these 
panels. These stars unambiguously determine the age of 
this population to be as old as the MW globular clusters 
(10 - 13Gyrs). 

Additionally information on the average age and metal- 
licity come from th e colour and magnitude of the core he- 
lium burning stars. iGirardi fc Salaria (|200ll ) show that the 
mean /-band magnitude of core helium burning stars is pri- 
marily dependent on age, whereas the mean [V — I) colour 



is primarily dependent on metallicity. While their study is 
focused towards understanding the evolution of red clump 
stars, the results they provide extend to old (low mass) he- 
lium burning stars as well. The top panels of Figure [lOl show 
the extinction-corrected V- and J-band luminosity functions 
of all stars in And II with r > 7.5 arcmins. The peak of stars 
in both panels is due to the HB. The dot-dashed lines corre- 
spond to the best fits of a Gaussian (which models the lumi- 
nosity profile of the HB ) and a power law (which models the 
slope of the RGB in this region). The peak magnitude of the 
Gaussians are the mean magnitudes of the HB stars. We find 
Vhb = 24.76±0.01 and I H b = 24.04±0.01, with correspond- 
ing dispersions of ay — 0.17 ± 0.01 and 07 = 0.14 ± 0.01. 
This implies My B = 0.69 ±0.06, Mf B = -0.03 ±0.06, and 
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Figure 11. Left panel: The metallicity distribution function (MDF) for stars in the top two magnitudes of the RGB with 10 > r > 
7.5 arcmins in Andromeda II (solid histogram). This radial range predominantly samples stars in the spatially extended component in 
Andromeda II. This was created by interpolating between 13 Gyr isochrones from VandenBerg et al. (2006) with BVRI colour- T e f j 
relations as described by VandenBerg & Clem (2003). The mean metallicity is [Fc/H] = —1.5 ±0.1 with frp e /jj] = 0.28 dex. Right panel: 
The same, but for RGB stars with r < 2 arcmins in Andromeda II (solid histogram), and interpolating between 9 Gyr isochrones. This 
radial range predominantly samples stars in the exponential component in Andromeda II. The mean metallicity is [Fe/H] = — 1.2±0.1 with 
°"[Fe/H] = 0.40 dex. The dashed histograms in both panels show the relevant MDF derived for the foreground population (r > 14 arcmins), 
scaled by area. 
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Figure 12. MDFs for the spatially extended and exponential components in And II (solid histograms; left and right panels, respectively). 
The dotted histograms show the MDFs in the outer and inner regions of And II for comparison. See text for details. 



(V — I) a — 0.72 ±0.02. The uncertainties in the magnitudes 
are dominated by the uncertainty in the distance modulus 
of And II. 

The relationships derived by iGirardi fc Salarisl (|200ll ) 
between the colour, magnitude, age and metallicity of core 
helium burning stars are shown in their Figure 1. By com- 
paring this with our observed values for the luminosity and 
colour of the HB stars we find that either these stars be- 
long to a very young, comparatively metal-rich population 
(< 2 Gyrs, [Fe/H] > —0.5), or they belong to a very old, 
metal-poor population (~ 13 Gyrs, —1.7 < [Fe/H] < —1.3). 



However, the blue HB stars show that this population must 
be old, favouring the more metal-poor estimate. 

A final constraint can be placed on the metallicity of the 
stars in this radial range by constructing a MDF from inter- 
polation of the position of stars on the RGB. We use exactly 
the same technique as before, adopting 13 Gyr isochrones for 
the interpolation procedure to ensure consistency with the 
HB age information. The result is shown in the left panel 
of Figure [TT] (for stars with 7.5 < r < 10 arcmins). Also 
shown in this panel as a dashed histogram is the MDF for 
stars with r > 14 arcmins, which is dominated by foreground 



14 McConnachie, Arimoto & Irwin 



stars. Clearly, the influence of foreground contamination on 
this MDF is minimal. The implied median metallicity of the 
(foreground corrected) MDF is [Fe/H] = —1.5 ± 0.1, with a 
dispersion of <T[f c /h] = 0.28 dex. The mean metallicity is the 
same as we derived in Section 3.1.1, although the MDF for 
stars at large radii is notably narrower than the overall MDF 
for And II (Figure Q. The excellent agreement between the 
age and metallicity estimates from the HB and the metallic- 
ity estimate of the RGB demonstrates good consistency in 
our analysis. 



3.5.2 Populations at small radii 

The bottom panels of Figure|9]show the CMDs for the inner 
few arcminutes of And II, which we expect is dominated by 
the exponential stellar component. However, approximately 
one-third of the stars are expected to belong to the spatially 
extended component. The RGB in this radial range is ex- 
tremely broad, and the morphology of the HB/red clump 
is very different to any of the CMDs examined previously. 
Some stars occupy the locus of the blue HB defined earlier: 
we count a total of 56 ± 8 such stars. However, our earlier 
analysis of the radial profiles shows that the expected num- 
ber density of HB stars at this radii is ~ 10 stars arcmin -2 
(Figure [8|, or a total of 100 HB stars in the area probed 
by these CMDs. Of these, ~ 43 % are expected to be blue 
HB stars. We concluded earlier that the HB stars traced 
an extended component; the expected contribution of blue 
HB stars from the extended component at r < 2 arcmins is 
therefore ~ 43 ± 7, which is consistent with the number we 
actually observe. We cannot rule out the possibility that the 
exponential component has a very weak blue HB, but this 
cannot contribute more than ~ 10 % of the total number of 
blue HB stars in And II. 

The lower panels of Figure [10] show the extinction- 
corrected luminosity functions for stars with r < 2 arcmins. 
The dashed lines show the best-fit Gaussians to the red 
clump/HB peaks. The mean magnitudes of the core helium 
burning locus are VheB = 24.73 ± 0.02 and IneB = 23.81 ± 
0.02, with dispersions of a v = 0.37 ± 0.05 and 07 = 0.34 ± 
0.02. These correspond to M$ eB = 0.66 ± 0.06, Mf eB = 
-0.26 ± 0.06 and (V - I) o = 0.92 ± 0.03. By compa r ing th e 
(V — I) a colour with Figure 1 of iGirardi fc Salarisl (2001), 
we imply a mean metallicity of [Fe/H] ~ —0.7 and a mean 
stellar age of ~ 9 Gyrs, with a plausible range given the un- 
certainties of 7 — 10.5 Gyrs. The mean /-band magnitude 
of the red clump is notably brighter than the mean HB lu- 
minosity measured previously, implying a stellar population 
that is significantly younger by at least 3 or so Gy r s. Thi s 
finding is in good agreement with |Pa Costa et alj {2000), 
who implied that there was a younger population in And II 
and that its age was probably in the range ~ 6 — 9 Gyrs. 

We derive an MDF for the inner regions in the same way 
as before, except we now use isochrones which are 9 Gyrs 
old to agree with the red clump age estimate. The MDF 
is shown in the right panel of Figure 111! and the dashed 
histogram shows the corresponding MDF for the foreground 
population scaled by area. The median metallicity is [Fe/H] 
= —1.2 ± 0.1, with a very broad dispersion of <7[ Fc / H ] = 
0.4 dex. Clea rly, the broad metallic ity d ispersion for And I I 
measured bv lDa Costa et all (|2000T ) and lCote etafl <|l999bl ) 



reflects the large dispersion present in the central regions of 
this galaxy, and which is not found at large radius. 

3. 5. 3 MDFs of the two components 

In our two component picture of And II, the MDFs shown 
in the left and right panels of Figure 11 are expected to pre- 
dominantly sample the spatially extended and exponential 
components, respectively. However, neither will give a clean 
sample of that component. 

To derive the intrinsic MDFs for the extended and ex- 
ponential components, we note that any MDF which sam- 
ples some radial range of And II will be a weighted sum of 
the two intrinsic MDFs, [MDF] a and [MDF] B , such that 

MDF = x[MDF] A /(l + x) + [MDF] b /(1 + x) (5) 

where x is the ratio of the number of stars in component A 
to component B, and is found by integrating the radial pro- 
files of the two components in the appropriate radial ranges. 
If we produce two MDFs (using the same age assumption) 
for two different radial ranges in And II, then we can solve 
these equations to find [MDF] a and [MDF]b (for that age 
assumption) . 

Figure 12 shows the intrinsic MDFs for the spatially 
extended component (left panel) and the exponential com- 
ponent (right panel) using this technique (under the age 
assumptions derived earlier). The dashed histograms show 
the MDFs for the outer and inner regions calculated previ- 
ously. The MDF for the extended component is unsurpris- 
ingly very similar to the MDF for the outer regions, given 
the complete dominance of this component at large radius. 
On the other hand, the MDF for the exponential component 
is significantly different to the MDF for the inner regions 
given the non-negligible contribution of the extended com- 
ponent in this radial range. The MDF for the exponential 
component implies a mean metallicity of [Fe/H] ~ —1, in 
reasonable agreement with the estimate from the colour of 
the red clump. Despite being significantly metal-rich, this 
component nevertheless appears to possess a large tail to 
very low metallicity. Interpretation of this feature is difficult, 
however, since age-metallicity degeneracies will undoubtedly 
have an effect, and the age of this component means that 
there could be contributions from Asymptotic Giant Branch 
stars masquerading as bluer RGB stars. 

We conclude that, in the two component model of 
And II, the spatially extended component is old (~ 13 Gyrs), 
has a mean metallicity of [Fe/H] ~ —1.5, and a dispersion 
of <T[ Fc / H ] = 0.28 dex. The exponential component, on the 
other hand, is younger (~ 7 — 10 Gyrs) and more metal-rich 
([Fe/H] ~ —1) with what appears to be a significant tail to 
low metallicities. 



4 DISCUSSION 

4.1 Two components versus smoothly varying 
gradients 

And II is a complex dwarf galaxy, both in terms of its struc- 
ture and its SFH. The inner region of the galaxy possesses 
very different stellar populations to the outer region. Ra- 
dial gradients in the stellar populations of dwarf galaxies 



The structure of Andromeda II 15 



ANDROMEDA II 



Old 

(extended) 



Intermediate age 
(concentrated) 



r (arcmins) 
To (PC) 
% light 
My 

L v (L Q x 10 6 ) 
M0,V (mags arcsec - 2 ) 
t (Gyrs) 
[Fe/H] 
C[Fe/H] 



0.3 ±0.1 
10.0 ±8.5 
1897 ± 1612 
~ 74 

i o q+0.3 

~ 26.1 
~ 13 Gyrs 
-1.5±0.1 
0.28 dex 



1.0 ±0.3 
2.2 ±0.8 
417 ± 152 

~ 26 
-11.2 
2.5 



+0.3 
0.2 
+0.6 
0.5 
24.8 



• 7 



10 Gyrs 
~ 1.0 
0.40 dex 



Table 1. Summary of properties for the two stellar components 
in Andromeda II 



have been observed many times before, most notably by 
lHarbeck et al.l (|200ll ) . Radial gradients can be produced due 
to the presence of two or more dynamically distinct stellar 
components with different stellar populations. In this sce- 
nario, the gradient arises as a result of the changing contri- 
butions of the two components as a function of radius. Al- 
ternatively, gradients could arise within a single dynamical 
component. For example, a centrally concentrated, chem- 
ically homogeneous gas cloud will form stars everywhere; 
however, more stars will be formed in the higher density 
central regions, which will enrich the central parts of the 
cloud more than the outer parts. Subsequent generations of 
stars at the center of the cloud will therefore form out of 
more enriched gas than stars in the outer parts, and so will 
evolve to occupy different loci in the colour-magnitude dia- 
gram after a fixed time, thereby producing a gradient in the 
stellar populations. 

For the specific case of the population gradients in 
And II, we strongly favor the first scenario of multiple dy- 
namical components. Compelling evidence in favour of this 
interpretation comes from the radial profiles shown in Fig- 
ure [8] The radial distributions of the red RGB and HB pop- 
ulations are dramatically different, and effectively possess 
different functional forms; the former is an exponential pro- 
file, whereas the latter is constant density out to large radius. 
It is very unlikely that two populations with such disparate 
spatial distributions could have formed out of the same gas 
cloud and still possess the same kinematics. If both popu- 
lations are in equilibrium with the dSph then they cannot 
both have the same kinematic properties while possessing 
such vastly different density profiles and still satisfy Jeans 
equation. If instead we argue that one or both of the pop- 
ulations is not in equilibrium, then we have no reason to 
expect the two populations will have the same kinematics. 

In addition to these considerations, if we adopt the as- 
sumption that And II consists of two distinct components 
then the overall radial profile of And II must equal the sum 
of the radial profiles of its two components. We find that a 
weighted sum of the HB and red RGB profiles fits the over- 
all radial profile of And II very well. We also find that the 
radius at which the central component begins to dominate 
is coincident with the radius at which MI6 noticed their pro- 
file for And II deviated significantly from a single component 
fit. In the Suprime-Cam data, the stellar populations inside 



this radius are very different to the populations at larger 
radii, and have different metallicities and mean ages. Given 
all of these factor, we conclude that the change in the stellar 
populations of And II as a function of radius is due to the 
presence of two, dynamically distinct, stellar components 
with distinct stellar populations. 

Table 1 presents a summary of the mai n properties we 
derive for the two components in And II. |Pa Costa et al.l 
(2000) first showed that the chemical abundance distribu- 
tion of And II could not be explained by a single component 
simple chemical enrichment model and instead required two 
contributing populations of different ages. The metal-poor 
population had a mean metallicity of log (< Z > /Zq) = 
— 1.6 and the metal-rich population had log (< Z > /Zq) = 
—0.95. The former outnumbered the latter by ~ 2.3 : 1. The 
metallicity estimates we derive for the two components in 
And I I agree well with the values implied bv lDa Costa et all 
(2000); in addition, we calculate the ratio of the metal-poor 
population to the metal-rich population to be ~ 2.8 : 1 (as- 
suming equal mass-to- light ratios), in good agreement with 
their model. We also find that the younger, more metal-rich 
population possesses a red clump; this is the first time this 
feature has been observed in a M31 dSph companion. The 
M31 dSphs are often characterised as having no significant 
intermediate age populations; we stress that these earlier 
results for And II and our own discovery of a red clump 
highlight that this characterisation is probably inadequate. 



4.2 Andromeda II and the Local Group dwarf 
spheroidals 

iTolstov et al.l ([2004.) were the first to show convincingly that 
the Sculptor dSph consisted of two, spatially, chemically and 
dynamically distinct populations. This is contrary to the 
prevailing treatment of dSph galaxies as single stellar com- 
ponent systems supported by velocity dispe rsion, embedded 
within a dark matter halo. Subsequently, Battagli a et al.l 
(2006) have shown that the Fornax dSph consi sts of three 
dynam ically distinct populations of stars, and llbata et al.l 
(2006) have shown the presence of similarly distinct pop- 
ulation s in Canes Venatic i. Howev er, equivalent stud ies of 
Leo I i|Koch et al.ll2007bl) , Leo II i|Koch et al.l l2007al ) and 
Carina (|Koch et all 120061 ) have not revealed any dynam- 
ically distinct stellar populations in these galaxies. In all 
cases where dynamically distinct components and/or radial 
gradients have been observed, the comparatively metal-rich 
component is shown to be more centrally concentrated than 
the metal-poor component. 

And II appears to bear some resemblance to Sculptor 
and Fornax. In all these galaxies, the more concentrated 
popul ations are you n ger th an the more extended popula- 
tions. iKawata et all (2006) examine the evolution of the 
baryonic component of dwarf galaxies in a cosmological sim- 
ulation and show that they can reproduce this general fea- 
ture as a result of multiple episodes of star formation: later 
episodes of star formation are spatially less extended than 
earlier episodes due to the dissipative collapse of the gas 
component. In this model, the populations form in situ, 
although this does not prevent the gas which formed the 
younger popula tions from being recen tly accreted. 

In contrast. [Battaglia et al.l l|2006l ) present tentative ev- 
idence for non-equilibrium kinematics in the most extended 
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(oldest) stellar component in Fornax. This might suggest 
that its formation was driven by mergers, perhaps by ac- 
creting an old dwarf galaxy or even a globular cluster which 
was subsequently de s troyed in the tidal field of F ornax (see 
IColeman et all 12004 IColeman fc Da Costal 120051 for tenta- 
tive evidence of past mergers in Fornax). Ilbata et al.l (|2006l ) 
also tentatively propose that the two components in Canes 
Venatici are not in equilibrium with the same dark matter 
halo. 

Other than a n overall velocity dispersion for And II 
based on 7 stars (ICote et al.l Il999ah . high quality global 
kinematics for And II are lacking. It is tempting, however, 
to speculate on the origin of the bizarre density profile of 
the extended component of And II. This component is ef- 
fectively constant density out to large radius (~ 10 arcmins 
~ 1.9 kpc) before truncating. This has not been observed in 
any other dwarf galaxy to date and is in contrast to the usual 
King, exponential, Plummer etc profiles of dwarf galaxies. 
Interestingly, numerical simulations of the tidal disruption 
of dwarf galaxies in MW-like potentials show some evidence 
of producing debris w hich is distributed rough ly evenly with 
radius (eg. Figure 3 of iPenarrubia et alJl200a ). Scaled down 
to dwarf galaxy potentials, the same physics may be respon- 
sible here; that is, the stars in the extended component of 
And II may be the debris from the accretion and eventual de- 
struction of an old stellar system in the potential of And II. 
Thus perhaps this component is a relic of the hierarchical 
formation of And II. 

While intriguing, detailed simulations of the tidal de- 
struction of stellar systems in the potential of a dwarf galaxy 
do not yet exist to study the above proposal in any de- 
tail. Indeed, a strong argument against this scenario is that 
the extended component of And II is the oldest compo- 
nent and contributes the majority of the light. As such, we 
would expect it to be the main body of And II and not 
accreted debris. The exponential component formed later, 
perhaps through a mechanism similar to that proposed by 
iKawata et al.l ([2009) . While this does seem the most natural 
explanation, we are forced to address the question of why 
the original body of And II has such an unusual and unique 
density profile. This is a fascinating question, and one which 
we are currently unable to answer. 



nents have very different stellar populations; the exponential 
component has an average age of ~ 9 Gyrs old, is relatively 
metal-rich ([Fe/H] ~ — 1) with a significant tail to low metal- 
licity, and possesses a red clump. The extended component, 
on the other hand, is ancient, metal-poor ([Fe/H] ~ —1.5) 
with a narrower dispersion o"[F e /H] — 0.28, and has a well 
developed blue HB. 

The extended component contains approximately three- 
quarters of the light of And II. Its unusual surface bright- 
ness profile is unique in Local Group dwarf galaxies, and 
implies that its formation and/or evolution may have been 
quite different to most other dSphs. The two component 
structure, however, is also observed in Fornax and Sculptor, 
and it may be that the centrally concentrat ed exponential 
compo nent formed in the way envisioned by IKawata et al] 
(2006) through the dissipative collapse of gas after the main 
component had formed. The chemo-dynamical structure of 
And II is clearly very complex and warrants spectroscopic 
studies of its metallicity and kinematic properties. It lends 
yet further support to the accumulating body of evidence 
which suggests that the evolutionary histories of faint dSph 
galaxies can be as every bit as complicated as their brighter 
and more massive counterparts. 
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5 SUMMARY 

We have used Subaru Suprime-Cam to obtain multi-colour 
imaging of the M31 dSph And II to a depth equivalent 
to earlier HST-WFPC2 studies but over an area 100 times 
larger. We identify various stellar populations, including a 
red clump. This is the first time this feature has been de- 
tected in a M31 dSph, which are often characterized as hav- 
ing no significant intermediate age populations. In the case 
of And II, this characterization is incorrect. 

We construct radial profiles for the various stellar popu- 
lations and show that the HB has a nearly constant density 
spatial distribution out to large radius, whereas the red- 
dest RGB stars are centrally concentrated in an exponential 
profile. We argue that these populations trace two distinct 
structural components in And II, and show that this as- 
sumption provides a good match to the overall radial profile 
of this galaxy. We go on to demonstrate that the two compo- 
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